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SUMMARY 



Condenser tubes with enhanced surfaces show promise for 
use in advanced marine condenser designs. In order to ade- 
quately predict the performance of enhanced tubes when operating 
in a condenser-bundle environment , considerable improvement in 
analytical method is necessary beyond the standard approaches 
recommended for use with conventional designs. 

This report describes a computer code of moderate complexity 
that is suitable for predicting the relative merits of the use 
of enhanced tubes in radial-flow condensers. The code accounts 
for vapor shear effects, condensate inundation, and noncondensi- 
ble gases, and predicts the local steam-side conditions at inter- 
mediate locations within the condenser bundle. 

Several computer experiments are described in which the 
impacts of enhanced-tube designs on overall condenser performance 
are investigated. The results serve to illustrate the importance 
of internal detail in such analyses and, in addition, they show 
that judicious use of enhanced tubes can lead to considerable 
improvements in condenser performance as well as the realization 
of heretofore unattainable design options. 
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PERFORMANCE RATING OF ENHANCED 



MARINE CONDENSERS 



INTRODUCTION 



Since the early days of naval steam propulsion, the sur- 
face condenser has evolved into a very reliable component of 
the steam propulsion machinery plant. This reliability has 
been achieved by providing generous design margins to insure 
thermal performance at full power. The penalties for this 
overdesign, however, are additional weight and volume (which 
must be carried around for the life of the ship) and crowded 
machinery rooms with poor accessibility. Although surface 
vessels and submarines have different design constraints, both 
types of vessels can benefit from a more compact condenser de- 
sign through improved heat -transfer. The actual dimensions 
of the condenser may have an impact on vessel size and dis- 
placement, the location of the propulsion machinery aboard the 
vessel, and even on vessel performance and cost. In addition, 
the future development of compact steam systems, both for main 
propulsion and secondary heat recovery purposes, will require 
the application of advanced technology to all system compo- 
nents, including the condenser. 

At the Naval Postgraduate School, the relative benefits 
of such advanced condenser designs have been analyzed using 
a modified, radial flow computer code (MORCON) to model the 
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heat transfer and fluid flow phenomena which occur within 
this type of heat exchanger. Refinements to this code and 
improvements in the models used to predict two phase heat 
transfer and pressure drop are continuing. 

Based upon the work accomplished to date, it is clear 
that present day smooth-tube steam condensers, when operating 
under typical conditions, are limited in thermal performance 
due to a large thermal resistance which occurs on the sea- 
water side (tube side) of the condenser. This resistance, on 
the average, is larger than the thermal resistances which oc- 
cur on the steam side and in the tube wall, and those due to 
seawater fouling and noncondensible gas buildup in the vicinity 
of the tube. Consequently, in present-day conventional con- 
denser designs, what happens on the steam side does not have 
to be modeled too accurately to arrive at reasonable designs.^" 
When using enhanced tubes, however, the enhancement on the sea- 
water side can increase the inside heat transfer coefficient by 
100 to 200 percent over the smooth-tube case. The outside heat 
transfer coefficient, on the other hand, is increased by only 
10 to 50 percent. In this situation, the thermal resistances 
on the inside and outside of the tube can be approximately 
equal. As a consequence, it becomes essential to understand 
and accurately model the detailed mechanisms which occur when 
steam condenses on these enhanced surfaces. 

‘'"It should be noted that, even with conventional designs, there 
can be locations within a condenser bundle where the steam side 
resistance is largest due to low vapor velocities and non- 
condensible gas pockets. 
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Heat transfer on the steam side of Naval condenser's 



often occurs in the presence of high vapor velocities and 
with severe condensate inundation from neighboring tubes in 
the bundle. It is also important to predict the resulting 
two phase pressure drop in order to trace the steam flow 
path through the condenser bundle and minimize losses in the 
saturation pressure and temperature across the bundle. Since 
comprehensive test data on this subject are virtually non- 
existent, and since advanced computer codes will depend 
strongly upon good steam-side models, additional experimental 
work is necessary. 

In spite of the tentative nature of existing theoretical 
steam-side heat transfer models, it has proved to be useful 
to conduct computer-based condenser performance calculations 
in order to provide comparative information regarding the rela- 
tive costs and benefits of enhanced- tube condensers. In addi- 
tion, the present theoretical basis is adequate to permit the 
identification of condenser tube-bundle regions that should be 
particularly sensitive to the introduction of internal and/or 
external enhancement. 

This report documents these computer-based calculations. 
Since steam-side conditions are especially important in such 
calculations, some preliminary discussion is devoted to this 
subject. This is followed by a brief general description of 
computer modelling of surface condensers. The computer code 
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used at NPS is then described and several case-studies are 
provided to illustrate the use of the code and to give some 
insights concerning the impact of enhanced- tube technology 
on condenser design. 

EFFECTS OF LOCAL STEAM-SIDE CONDITIONS 

Single- tube theories and experiments have shown, not 
surprisingly, that the thermal performance of a condensing 
tube may be dramatically affected by the thermo-fluid state 
of the surrounding two-phase multi-component flow. Thus the 
condensate loading and velocity of the vapor, its temperature, 
and its non-condensible gas content all have important influ- 
ences on the thermal resistance on the steam side of the tube. 
These influences lead to further complications in the analy- 
sis of condenser tubes performing in a bundle environment and 
the overall bundle performance is strongly coupled to the 
bundle inlet conditions, individual tube performance, and 
the general nature of the tube distribution within the bundle. 
Several of the main factors that are likely to be significant 
in bundle performance are discussed in the following para- 
graphs. 

LOCAL STEAMnSIDE PRESSURE VARIATIONS 

Frictional pressure losses on the steam side of a con- 
denser bundle are detrimental from a system point of view 
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because these losses are directly chargeable against the 
thermal efficiency of the power system. In addition, the 
vapor temperature, which is a major factor in the heat trans- 
fer process that accompanies condensation, is strongly de- 
pendent upon the pressure in the vapor. This dependency is 
accentuated at low condensing pressures, varying from about 
7°F/in. Hg at 6 in. Hg condensing pressure to 12°F/in. Hg at 
3 in. Hg. (In other words, at about 3 psia a loss of 0.1 psia 
in pressure will lead to a 1.3°F loss in vapor temperature and 
at 1.5 psia the loss would be about 2.4°F.) Consequently, 
when condenser tubes are closely packed in a bundle the re- 
sulting high vapor velocities and correspondingly high pres- 
sure losses can severely limit the bundle performance. 

VAPOR VELOCITY EFFECTS ON STEAM-SIDE THERMAL RESISTANCE 

High vapor velocities can have at least one effect that 
is beneficial. This is the shearing force that the vapor ex- 
erts on the condensate film which, at relatively high vapor 
velocities, can be significant. A number of investigators 
Cl, 2, 3, 4] have contributed to the understanding of the 
effects of vapor shear. In particular, Fujii, et al. [1] 
have recently proposed experimental correlations of the form 
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where Nu is the mean Nusselt number, Re T is the "two-phase" 
Reynolds number (based on vapor velocity, tube outside di- 
ameter, and kinematic viscosity of the condensate ) , and Nu q 
is the standard Nusselt number for the ideal (zero shear) 
case. The empirical constants and a are in the ranges 

1.13 < C 1 < 1.24 

.195 < a < .2 

depending upon whether the tube thermal conditions are most 
appropriately described by constant temperature or constant 
heat flux conditions. The larger values are for the constant 
heat flux case. According to Fujii, et al. , Eq. (1) is valid 
in the range 3.3 > Re_ /Nu > 0.28. For smaller values of this 
parameter, they recommend the use of a slightly reduced value 
of the standard Nusselt number that is given by 

Nu = 0.96 Nu 
IU o 

In a typical naval condenser, vapor velocities well above 
100 ft/s are not uncommon and two-phase Reynolds numbers on 

g 

the order of 10 are typical. For a nominal Nusselt number 
value of Nu q = 400, the resulting vapor shear effect is given 
by Nu m /Nu Q ~ 1.5. Thus, the vapor shear effect is clearly 
important in the range of operation of naval condensers, and, 
when enhancements of various kinds are applied, this effect 
may become a dominant factor in overall condenser performance 
predictions . 
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Correlations of the form of Eq. (1) appear to be ade- 
quate for the prediction of shell-side film coefficients in 
horizontal and downward vapor flow without condensate inun- 
dation. For upward vapor flow, however, the effects of vapor 
shear on film thickness appear to be much larger, and at very 
high velocities the shell-side film coefficient, normally low 
due to poor condensate removal, can be significantly increased 
- even to the extent that the influence of gravity is overcome 
by vapor shear effects. 

From the point of view of enhancement, simple reasoning 
would lead to the expectation that improvements in heat trans- 
fer due to vapor shear will be less significant at high ex- 
ternal enhancement conditions. This follows from the argument 
that enhanced tubes will carry thinner faster-moving conden- 
sate films (other conditions being equal) that are less 
sensitive to shearing effects acting at the vapor interface. 
Equation (1) reflects this effect in- that (since a < .25 in 
all known cases) the ratio Nu /Nu varies inversely with Nu . 
The vapor shear effect might also be reduced by some of the 
previously-mentioned external enhancement schemes employing 
extended surfaces that would trap or otherwise "protect" the 
condensate. With high levels of internal enhancement, on the 
other hand, failure to provide for adequate condensate removal 
will be somewhat compensated by the shearing of the relatively 
thick condensate film at high vapor velocities. 
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Other vapor velocity effects, such as rippling, induce- 
ment of turbulence and stripping of the condensate film, are 
relatively poorly understood. A particularly important area 
in need of further research is the consideration of the en- 
trainment of condensate droplets in fast-moving vapor flows. 

As will be discussed in the subsequent section, condensate 
droplets as they traverse the condenser can have a profound 
effect upon the performance of downstream tubes. To the ex- 
tent that the paths of these droplets and the vapor flow are 
coupled through viscosity, an additional dependency upon vapor 
velocity must be considered. 

CONDENSATE INUNDATION 

The true behavior of the condensate generated within a 
tube bundle is far from the ideal model of vertically-dropping 
laminar sheets hypothesized by Nusselt. In fact, the conden- 
sate in most condenser environments will be in many forms 
(droplets, rivulets, sheets, sprays) and will flow in many 
directions. It is highly unlikely that any quantity of con- 
densate will traverse a condenser without losing its identity 
by impinging upon other tubes in the bundle. The effect of 
this impingement is referred to here as condensate "inunda- 
tion" or "rain" and it is probable that a true theoretical 
model of the phenomenon would of necessity be stochastic in 
nature. In the absence of such analytical sophistication, a 
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number of empirical models have been developed that are useful 
when applied within their range of credibility. 

Strict application of the classical Nusselt theory, to- 
gether with its restrictive assumption of a laminar film fall- 
ing under the influence of gravity alone, leads to the conclu- 
sion that the average condensate film heat transfer coefficient 
for a vertical row of n tubes, h , is given by [5] 

\ = n -0.25 (2) 



where h is the single-tube film coefficient. Numerous inves- 
o 

tigators have compared this simple expression with experimental 
results [6, 7, 8, 9, 10, 11] and it has become clear that in 
reality the results obtained from Eq. (2) reflect an over- 
prediction of the detrimental effects of inundation. In fact, 
a modified form of Eq. (2) may be written as 




h 



o 



n 



-s 



(3) 



where empirical evaluation of s has led to values in the range 
.07 < s < .25 [7], Within this range, for example, a column of 
30 tubes may be judged to have an average film coefficient that 
is 79% to 43% of the single- tube value. It should be noted 
further that the actual film coefficient for the n-th tube in 
a column, h , may be found using Eq. (3) with the result that 
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n _ ( 1-s) / , x ( 1-s) 

— = n - (n-I) 



(4) 



The range of s leads to a range in h n of from 75% to 32% of 
h Q : clearly the range of uncertainty in s is unacceptable. 

Among the chief discrepancies in the assumed constancy 
in s in Eqs. (3) and (4) are: 

1. The effect of inundation is definitely dependent 
upon the vapor velocity and direction. 

2. Bundle tube layout (tube spacing and pattern) will 
influence condensate inundation. 

Thus, even if the functional relationship given in Eq. (3) is 
authentic, the exponent s must show the dependency 

s = fct (steam velocity vector, bundle geometry) . 

The first of these is indicated nicely by Brickell [ll] in 
which a general decrease in s with steam velocity is shown. 
These results are reproduced here in Fig. 1. It is also ap- 
parent that at low steam velocities the exponent s is larger 
(more inundation effect) for both upward and downward steam 
flow than it is for horizontal steamflow. This directional 
effect is less noticeable at high steam velocities. The in- 
fluence of tube separation upon inundation effects is not 
clearly discernable from available data, largely because of 
the coupling of vapor shear effects with those due to inunda- 
tion. It is possible to speculate with some confidence, how- 
ever, that closely-spaced tubes are more influenced by 
inundation than are the tubes in relatively disperse bundles. 
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One important influence of tube spacing in low velocity 
steam flow is that of side-drainage whereby condensate gen- 
erated from above may, due to surface tension effects, pro- 
ceeds laterally to adjacent tubes rather than directly down- 
ward. Eissenberg [7] has contributed a means of quantifying 
this effect with the result that Eq. (3) may be modified as 
follows: 

h 

jp = (1 - F d )n~ S + F d [0.6 + 1/2 (2n)~ S ] (5) 

o 

where F d depends upon tube spacing and approaches unity for 
densely- packed staggered tube banks and zero for disperse tube 
layouts. Equation (5) represents the state-of-the-art in the 
prediction of condensate inundation effects and its effective- 
ness in doing so is, of course, highly dependent upon the 
selection of correct values for s and F d « 

The extent to which enhanced tubes are susceptible or 
resistant to inundation effects is a matter that is far from 
settled. (Indeed, the issue is unresolved as far as smooth 
tubes are concerned, when bundle effects are taken into ac- 
count.) There are some experimental results [12] that indi- 
cate that smooth tubes operating in the dropwise mode of 
condensation are not greatly affected by inundation. (This 
observation was made at very low vapor velocities.) Prelimi- 
nary indications are likewise encouraging for film condensation 
on enhanced tubes where single- tube experiments [13] have indi- 
cated values of s [Eq. (3)] in the range of 0.04 to 0.10. 
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NON-CONDENSIBLE GASES 



The effects of non-condensible gases can be classified 
into one of two categories: 

1. The introduction of an additional local thermal 
resistance due to the propagation of gases towards 
a condensing tube surface under the influence of a 
gradient in partial pressure, and 

2. The cumulative effect of gas blanketing where uneven 
rates of condensation in a condenser bundle eventu- 
ally lead to regions where tubes are inoperative in 
a condensing role. 

In naval condensers [14] as well as in shore-based power 
condensers [8], non-condensible gases contained in the inlet 
flow are usually on the order of 100-200 weight-parts-per- 
million. This ratio increases in the condenser, as steam is 
extracted, but even at a worst-case condition (20,000 wppra, 
say) the local overall heat transfer coefficient will be re- 
duced by only about 15% and this in a very small region of 
the condenser. Thus, although a number of theories are avail- 
able to account for the local effects of non-condensible gases 
[15] in naval condensers, there appears to be little need for 
extreme precision in the calculation. 

The effect of non-condensible gases on the external re- 
sistance of enhanced tubes is a matter of considerable uncer- 
tainty. Cunningham [16] has recently reported on single-tube 
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experiments in which fine-finned tubes were found to be less 
affected by non- condensible gases than were smooth tubes 
under corresponding conditions. In any event, it is ex- 
pected that because of the low concentrations of non-condensible 
gases under naval condenser conditions their effect will be 
small unless their presence in some way affects the enhancement 
process itself. 

The more- important aspect of non-condensible gas effects, 
crucial in all condenser applications, is the design of tube 
bundles such that non-condensible gases do not accumulate. 

In large tube bundles it is often necessary to incorporate 
special inlet baffling and steam lanes that allow the steam 
to penetrate deep into all parts of the tube nest in order to 
avoid gas blanketing. The compensation of condenser bundle 
designs for this purpose has often been experimental in na- 
ture, involving large-scale and extremely expensive tests, 
and is probably a major reason why condenser cross-sections 
do not seem to change drastically from generation to genera- 
tion. 

In recent years, however, the ascendancy of computational 
capabilities has led to the feasibility of the comprehensive 
computer codes necessary for the prediction and alleviation 
of gas blanketing. One such computer code has been developed 
in the UK (such codes are not, in our experience, available 
to the public) and has been used to great advantage in the de- 
sign and improvement of shore-based power condensers [8]„ 
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SUMMARY 



The performance of a condenser tube within a tube bundle 
is strongly influenced by the history of the vapor/gas mixture 
prior to arrival at the tube in question. Conversely, the 
performance of any given tube can have a profound effect upon 
the thermal behavior of tubes located downstream. The com- 
plexity of these interactions may lead to a degree of depres- 
sion in the mind of the researcher who would predict overall 
bundle performance. On the other hand, these inadequacies in 
understanding may be viewed as challenges and, even in our 
present state of ignorance, we can answer questions in a 
relative way. For instance, qualitative trends due to aug- 
mentation may be identified and relative improvements due to 
augmentation may be quantified within a precision of perhaps 
10 - 20 %. 

Just as a condenser is larger and more complex than a 
single tube, it is to be expected that the necessary compu- 
tational and experimental scales will be larger for condensers 
than they are for single tubes. To justify such investments, 
and to guide the efficient allocation of resources, it is 
necessary to use what tools we have at hand to investigate 
potential payoffs and pacing problems. The following sections 
describe some preliminary results obtained from computer 
models of condenser bundles. 
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COMPUTER MODELING 



Analytical models for the estimation of condenser per- 
formance may be classed as follows: 

1. Standard methods of the ''black-box'' type that are 
used for preliminary - and sometimes, it is sus- 
pected for final - design purposes. These methods 
include those recommended by the Heat Exchange In- 
stitute (HEI [17]) and British Electrical and Allied 
Manufacturers' Association (BEAMA [18]). 

2. Intermediate computer codes that, to varying extents, 
include the details of tube-to-tube flow within the 
condenser. For purposes of classification, these 
codes are restricted to the analysis of condenser 
designs in which the vapor/gas mixture may be pre- 
sumed to flow in a fixed direction (.such as radial, 
horizontal, upwards, or downwards). The code MORCON 
is of this class and is described in some detail be- 
low. 

3. Multi-dimensional codes utilizing finite element or 
finite difference methods for the solution of the 
governing conservation equations. Because codes in 
this category are proprietary (in all cases, to our 
knowledge) , little is known concerning the extent to 
which they exist and work. The CERL code mentioned 
above [8] appears to be quite successful in application 
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and there are some indications [19] that (J.S. manu- 
facturers possess a similar capability. 

The black-box methods are extremely limited to their 
application to all but the most standard condenser designs. 

In essence, only the overall heat transfer coefficient, U, is 
predicted by these methods and even this is subject to con- 
siderable uncertainty (HEI predictions of U have been ob- 
served to differ from measured values by from 100% on the 
high side to 50% on the low side [20]. These methods of pre- 
diction are based upon the assumption that the overall thermal 
resistance is dominated by the internal (tubeside) resistance 
- they are limited to cases of moderate cooling water velo- 
city - and the allowances for fouling are subject to consid- 
erable criticism on physical bases. No means are provided 
for accounting for the distribution of thermal resistances . 
across individual tube sections or between tubes in a bundle. 
It is clear that major advances in condenser designs cannot 
be based on these black-box methods. 

At the other extreme of the spectrum of sophistication 
lie the multi-dimensional computer-based methods. The basic 
computational difficulties in the development of these com- 
prehensive computer codes are largely ones of scale. It is 
quite difficult, yet most important to include two-dimensional 
capabilities that allow calculations in the plane of flow 
normal to the axis of the condenser. The third dimension, 
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steam/gas flow in the direction of the tube axis, may be ob- 
tained with reasonable accuracy by combining two-dimensional 
slices. Even so, the complexity of the two-dimensional cal- 
culation is sufficient unto itself in that it is often neces- 
sary to compute the two-phase multi-component flow in the 
passages between thousands of tubes. As a result, it is 
common to combine groups of tubes into "mini-bundles". Within 
these mini-bundles the equations of mass, momentum, and energy 
are satisfied by means of finite-difference expressions for 
the convective, conductive, and source effects. The source 
terms in these expressions require much of the approximation 
necessary in less-complicated codes in terms of interaction 
effects such as vapor shear and condensate inundation. Thus, 
the extension to multi-dimensional computations does not ob- 
viate the need for additional insights into the various bundle 
interactions mentioned in the previous section. The chief 
benefit of such extensions appears to be the ability to pre- 
dict the performance of poorly designed condenser bundles 
especially with respect to air-blanketing. 

Multi-dimensional computer codes can be extremely valuable 
for condenser performance predictions and may even be essential 
to the development of truly optimized condensers. Large re- 
source investments are required for the development of such 
codes and questions inevitably arise concerning their cost- 
effectiveness. This issue is not developed further here except 
to say that, because of it, our emphasis has been at the 
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intermediate level of analytical sophistication where codes 
such as MORCON are found. 

THE MORCON COMPUTER CODE 

MORCON (Modified Oak Ridge Condenser) is a modified ver- 
sion of the 0RC0N1 code developed at the Oak Ridge National 
Laboratories [2l], Some of the features and limitations of 
MORCON are summarized below: 

a. MORCON allows the calculation of radial flow in cir- 
cular bundles. (Other unmodified versions of 0RC0N1 
have been developed for downward and horizontal-flows.) 

b. Correction of the condensate film coefficient for the 
effects of rain, including side drainage, is accom- 
plished by the means previously described. As has 
been pointed out, the correction only accounts for 
the distribution of condensate rain upon lower tubes 
in a gravity-dominated flow and does not account for 
the interaction of vapor shear effects with those 
due to inundation. 

c. Non-condensible gases (up to four species) are taken 
into account by the use of the Colburn analogy [22] 
together with an empirical correlation credited to 
Eissenberg in [21], 

d. MORCON allows for the effects of heat transfer en- 
hancement of the inner and/or outer film coefficients. 
The influence of enhancement on tubeside and shell-side 
pressure drops is also included. 
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e. MORCON, as in all other known cases, includes foul- 
ing effects only as a constant resistance term. 

f. All tube-side and shell-side fluid properties are 
continuously calculated as they change with pres- 
sure and temperature throughout the condenser. The 
latter quantities vary due to inlet, outlet, and 
tube- to- tube flow cross-section variations, as well 
as viscous losses. 

g. Bundles can be segmented into as many as six compu- 
tational sectors. Steam flow entering each sector 
is adjusted until the shell-side pressures at the 
exits of each sector (which is also the air cooler 
inlet) are matched within a specified limit. The 
sector arrangement for the circular bundle with 
radial flow is illustrated in Fig. 2. 

h. The overall heat transfer coefficient for each row 
in each sector is calculated on the basis of an 
"average" tube. This is a basic limitation of the 
code that is somewhat relieved by the ability to 
calculate performance by sectors. 

i. Vapor shear effects are included in the manner pre- 
viously described. To the extent that the Fujii 
correlation does not include the interaction of 
vapor velocity with condensate inundation, the re- 
sults are suspect for other than gravity-dominated 
(low vapor velocity) flow. 
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air removal point 




Figure 2 . Sketches showing 
sector arrangements and 
quantities necessary for 
inundation calculations 
[ 21 ]. 
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DESCRIPTIVE SUMMARY 



In this section the organization and flow of MORCON are 
described. A listing of the code is provided in Appendix A. 
Because MORCON is in large part identical to the ORCONl code, 
the treatment here will be brief with details confined to the 
few modifications made to ORCON 1. For a more complete des- 
cription of the basic code, the reader is referred to [21]. The 
subroutines contained in the code are listed below: 

MAIN (not labeled) 

SECALC 

COOLEX 

HETTRN 

ADJUST 

PRSRDP 

DFSVTY 

INPUT 

0UTP1 

ZERO 

Subroutine DFSVTY computes the mutual diffusivity of the vapor/ 
gas mixture. Various function subroutines are included to cal- 
culate specific heats, thermal conductivities, viscosities, and 
densities of the gas, vapor, and condensate as required. Func- 
tion subroutines are also provided to calculate vapor saturation 
pressures and temperatures and the latent heat of vaporization. 
The purposes of subroutines INPUT, 0UTP1, and ZERO are indicated 
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by their names and these will receive no further discussion here. 
The subroutine calling protocol is indicated in Table I, below. 

Table I. Subroutine Calling Protocol 



Calling 

Sub- 

Routine 


Called Subroutine 


SECALC 


COOLEX 


HETTRN 


ADJUST 


DFSVTY 


PRSDRP 


INPUT 


0UTP1 


MAIN 


X 


X 




X 






X 


X 


SECALC 






X 






X 






COOLEX 






X 






X 






HETTRN 










X 









Subroutine MAIN 

Subroutine MAIN is the main control program. Input values 
are taken into this subroutine and the tube bundle cross-sectional 
geometry is determined. This geometry is defined by the specified 
number of tubes, their outside diameter, and the spacing between 
them. The tubes are arranged radially such that the outer row of 
the bundle is complete. The tube array is staggered with the 
radial spacing between rows equal to (/3/2) times the specified 
circumferential spacing between tubes. The void diameter is de- 
termined by the space necessary to accommodate the number of tubes 
that are insufficient to complete an outer row. With the tube 
sheet entirely filled (zero void diameter) , the specification of 
additional tubes will result in the removal of a number of tubes 
from the void sufficient to complete an additional outer row. 
(Because of the difference in performance of inner and outer-row 
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tubes, this method of allocating tube cites leads to highly 
non-linear and sometimes confusing results when tube count is 
used as an independent design variable.) 

Bundle symmetry is assumed about the vertical centerline 
of the circular bundle and the remaining semi-circular tube 
array is divided into six sectors. Each partial row in each 
sector is defined by the central tube in that string and the 
program calculates the number of tubes in a vertical row above 
this "mean” tube. This is needed to provide for inundation 
calculations. 

Following these preliminary arrangements, subroutine 
SECALC is called. The return information is the condition of 
the vapor/gas mixture as it leaves the main condenser bundle 
and enters the air cooler. The performance of the cooler is 
then computed by means of a call to subroutine COOLEX. 

The final overall vapor/gas conditions, as returned from 
COOLEX, are transmitted to subroutine ADJUST. If the exit 
fraction of the current design is equal to the input specifi- 
cation (within five percent error) the return from ADJUST is 
followed by final output actions. If the specified exit frac- 
tion is not obtained, the code will adjust either the inlet 
steam flow or the tube (bundle) length, as selected by the 
user. The updated design is then returned to SECALC for evalu- 
ation and the above procedure is repeated until convergence is 
obtained or the iteration limit (set to 10) is exceeded. 

Figure 3, below, summarizes the function of subroutine MAIN. 
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Figure 3. Simplified flow diagram of main calling subroutine. 
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Subroutine SECALC 



Subroutine SECALC, with calls to subroutine HETTRN, de- 
termines the sector-by-sector and row-by-row performance of 
the condenser. Calculation is begun in the first (uppermost) 
sector at the first (outermost) row. Using the central tube 
of this and subsequent row segments, subroutine HETTRN is 
called to determine the thermal performance of the current 
"man" tube. With information returned from HETTRN, subroutine 
PRSDRP is called to determine the pressure of the vapor/gas 
mixture entering the next row and the steam quantity and con- 
ditions for this next row are returned to HETTRN. 

The outer loop in SECALC covers the six sectors while the 
inner loop covers the rows (the number of rows is determined 
by the input geometry) . Following the completion of both 
loops, the steam-side pressures at the outlet of each sector 
are each compared with the mean of the outlet pressures. The 
total inlet steam flow is then adjusted so that sectors with 
excessive steam-side pressure drop receive less flow and con- 
versely for those sectors with relatively higher outlet pres- 
sures. The calculation is then repeated in its entirety (all 
sectors, all rows) until the bundle outlet pressure is uniform 
to within one percent. 

Figure 4, below, illustrates the function of subroutine 
SECALC. 
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Figure 4. Simplified flow diagram for subroutine SECALC. 
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Subroutine COOLEX 



Subroutine COOLEX is a simplified version of subroutine 
SECALC in which calls to HETTRN are used to obtain row-by-row 
performance of the cooler tubes. No steam-side pressure 
balance (such as that performed in SECALC) is attempted in 
COOLEX. 

The user specifies the percentage of the total number of 
condensing tubes to be used in the cooler section. These 
tubes are allocated to a vertical staggered array with cross- 
flow dimension initially equal to the bundle void diameter. 

If this tube cite allocation results in crowding of the cooler 
tubes (as indicated by vapor/gas mixture velocities in excess 
of 150 ft/s) the cooler is widened (and proportionately 
shortened) as necessary. 

Subroutine HETTRN 

This program contains the basic heat transfer algorithms 
and, as noted above, is called by SECALC for each mean tube 
in each row of each sector. Following a summary of the pro- 
cedures followed in HETTRN, a few of the more- important heat 
transfer algorithms are discussed in some detail below. 

Because the thermodynamic and physical properties of the 
cooling water are allowed to depend upon temperature in HETTRN, 
a basic iteration is required to compute the overall heat 
transfer coefficient. This iteration begins with an estimate 
from which, with local steam conditions and cooling water inlet 
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temperature and flow rate, the cooling water outlet tempera- 
ture is calculated. Cooling water properties are then calcu- 
lated based upon the mean of the inlet and outlet temperatures 
of the cooling water. 

Following the determination of cooling water properties, 
HETTRN then computes, in order, the internal thermal resistance 
(based upon the Di ttus-Boelter correlation with internal en- 
hancement taken into account) and the external condensate 
film thermal resistance based upon the Nusselt correlation 
(see below) . The external resistance is then corrected for 
vapor shear and condensate rain effects. These resistances, 
together with tube conduction resistance and a constant 
fouling resistance (input by the user) , are then summed to 
develop an interim value of total resistance to heat transfer, 
exclusive of the resistance of the non-condensible gas film. 

The thermal resistance of the non-condensible gas film is 
calculated and various key variables for the current iteration 
are determined. Convergence of the basic loop is based upon 
agreement between current and previous values of film tem- 
perature drops. Agreement to within 0.01 degrees Fahrenheit 
is required in the case of the liquid film temperature drop 
and if non-condensible gases are present the iterations con- 
tinue until the calculation of the gas film temperature drop 
is steady to within 0.001 degrees Fahrenheit. Following con- 
vergence, values are transferred to SECALC and common storage 
as necessary for future use. 
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Condensate Film Thermal Resistance 
In MORCON the standard Nusselt relationship for the 
condensate film resistance is corrected for vapor shear ef- 
fects according to the correlations of Fujii et al. [l], On 
the assumption that the thermal resistance on the steam side 
is relatively small, the constant heat flux case is taken to 
be most appropriate. (Within experimental uncertainties the 
two cases - constant heat flux or uniform wall temperature - 
are not appreciably different.) Eq. (1) becomes 



Nu 

NT 2 - !- 24 R «L V 

o 



( 6 ) 



where Nu q is the Nusselt relationship given by 



Nu » 0.725 
o 



gp 2 d_ 3 h 






O fc g 



k lUl a T 



C7) 



Following correction of the mean Nusselt number for 
vapor shear effects, this value is further corrected for 
condensate inundation. Equation (5) is applied with a value 
of s = 0.2 so that the average film coefficient for a ver- 
tical row of m tubes is given by 



h _ o 2 

££ = 0.6F d + (1 - Q. 5647F d )m (8) 

This expression is evaluated for m = n and m = n-1, where 
n is the number of condensing tubes vertically above the 
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h„/h. 



tube in question, and the heat transfer coefficient for that 
tube is given by 

h = n h - (n-l) h « (9) 

n n n— x 

The selection of the exponent s = 0.2 in these calcula- 
tions is based upon the comparison with experiment shown in 
Fig. 5 and is also in consonance with the observations of 
others [7, 11, 23] (see also Fig 0 1) . 




Figure 5. Tubes with 
inundation and side 
drainage; comparison of 
theory with experiment. 
Data of Ferguson and 
Oakden ("Heat Transfer 
Coefficients for Water 
and Steam in a Surface 
Condenser," Trans. Chem. 
Eng. Congress (.World 
Power Conf.), v.3; 640 
(1936)) as reported in 
[7]. 



Non-Condensible Gas Film Coefficient 
The process of heat transfer through the non-condensible 
portion of the gas/vapor mixture is modelled as one due to 
mass transfer of the vapor through the gas. Sensible heat 
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transfer across the gas film (usually about 0.1% or less) is 
neglected. A mass transfer coefficient, K^, is defined such 



that 



K, 



V T = 

h fg AP sf h fg 




( 10 ) 



where AP^ is the drop in vapor partial pressure in the 
film, AT is the corresponding temperature change, is the 

latent heat of vaporization, and h^ is the heat transfer co- 
efficient for the process across the non-condensible film. 

The drop in vapor partial pressure across the gas film 
is taken equal to the rise in gas partial pressure across 
the film so that 



AP - = P - P , 
sf gc gb 



( 11 ) 



where P^ c and P^ are the gas partial pressures at the con- 
densate film interface and the external mixture conditions, 
respectively. Combining Eqs . (10) and (11), the gas film 

temperature drop may be written: 

K„h, 



AT = — ^ 3 - P rrr* ~ P 



s c ‘ gM H 



( 12 ) 



The partial pressure P ^ is known at each point in the calcu- 
lation and the rate (9T/9P) is calculated from the saturated 
steam line. The condensate interface gas partial pressure 
is taken from the empirical relationship suggested by Eis- 
senberg and Noritake [24]: 

C . 

P = -1 
gc k g 



32 



where 



exp (.54 - .544 In Re) 



D = 



( s c > 



2/3 



with G the local mass velocity, Sc the local Schmidt number, 
and Re the local vapor Reynolds number. With this correla- 
tion, Eq. (13) becomes 



AT 



= c j h f g , _ P (*1 

h g F gb \ 3P 



(13) 



The gas film temperature drop may be further defined in terms 
of heat transfer coefficients by the following heat balance: 



h AT = U AT 
g g ra 



(14) 



where U g is the overall heat transfer coefficient belonging 
to the log-mean temperature difference AT m . That is, if U is 
the heat transfer coefficient without non-condensible gases 
(known, at this point in the calculation) then: 



I = A + I 

U h rr U 

g g 



(15) 



Combination of Eqs . (13) - (15) with the elimination of AT and 

h g leads to the following quadratic relationship for the over- 



all heat transfer coefficient: 



U 



g - 




1 + 



gb 

AT 

m 



(n) 



h- c .u 

a + = o 

9 AT m 



(16) 
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All the quantities necessary to determine the coefficients 
in this relationship are available (for the current iteration) 
so that a solution for U is forthcoming. Solutions for h and 

g g 

AT follow according to Eqs. (15) and (14) . 

Subroutine PRSRDP 

As the name implies, the purpose of this subroutine is 
to compute the pressure drop on the steam side between adjacent 
rows of tubes. In MORCON (as in 0RC0N1) the friction factor is 
given by [24] : 

f = 0.102 + 52.2 (17) 

Re 



where f is related to the pressure drop by 



or 



4f 



AP 



AP = 2f G 2 

P 



The application of the correlation of Eq. (17) is limited to 
moderate mass velocities CG < 1500 lb/hr - ft 2 ) and, as such, 
may lead to significant errors when used to predict steam-side 
pressure drops in some advanced marine condenser designs. 



SAMPLE RUN 



In this section a sample run is described in order to il- 
lustrate the necessary user procedure and to provide a base- 
line case for subsequent analysis of the effects of enhancement. 
The tubesheet layout for the baseline design is shown in Fig. 6. 
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Figure 6. Assembly and tubesheet layout for baseline 
condenser . 
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and general specifications for this design are listed in 
Table II. 

Table II. Design and Performance Specifications for 
the Baseline Condenser 

Inlet steam: 189,003 lb m /hr saturated 5 in.Hg (133. 7°F) 

33.7 lb /hr noncondensibles (air) 
m 

Tubes: 1646 tubes in each of two passes, 6% in cooler 

70-30 Cu Ni, 5/8 OD, 16 BWG (0.065 wall) 
Pitch/OD 1.35 , , 

Fouling Resistance 3.3x10 hr ft“ °F/Btu 

Cooling water: Inlet temperature 66.1°F 

Velocity 8 ft/s (7,900 gpm) 

Salt concentration 3.5% by weight 

User Supplied Data 

Subroutine INPUT reads, writes, and processes the data 
provided by the user. These data are supplied via seven 
cards, as illustrated in Fig. 7. The first five of these 
cards are identical to those described for 0RC0N1 [21] but 
the input nomenclature is repeated here so that the reader 
may proceed (it is hoped) without reference to that document. 
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Figure 7. Card deck layout for input data. 



Following is a list of definitions of the input para- 
meters. A listing of the input cards for the baseline case 



is given in Appendix B (compare Table II) . 
Card 1, (FORMAT: II, 19A4, A3) 



IR 


Indicator for cards to be read (column 1 
only) . If IR = 0, 4, or >5, then all seven 
cards will be read in proper order. 


INDENT 


All the remaining columns of the first card 
are available for identification or notes, 
79 columns 



Card 2, (F0RMAT : 7F10.0, 215) 



ANTP 


Total number of tubes in the bundle (includes 
cooler) . 


PRCCLR 


Percent of tubes in the cooler. 


SDD 


Ratio of tube circumferential center- to-center 
spacing to tube diameter 


SDDMIN 


Ratio of tube spacing to tube diameter for 
cooler section. 


0D 


Outside diameter of all tubes, in. 


XW 


Wall thickness of tubes, in. 


SKW 


Thermal conductivity of tube material, 
Btu/hr-f t 2 -°F/ft. 


NRWSP 


Number of rows of condenser coolant inlet tem- 
peratures to be read (use only to read punched 
card inlet coolant temperatures, set = 0 other' 
wise) . 
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NRWSCP Same as NRWSP except for cooler. 



Card 3, (F0RMAT: 8F10.3) 



WBI 


Total coolant flow to condenser (cooler in- 
cluded) , lb/hr. 


VELBIP 


Coolant velocity, ft/sec. 

(Note: WBI and VELBIP are exclusive variables; 
one of them should be set to zero. If not, 
then VELBIP will override WBI) . 


CBI 


Coolant salt concentration, fraction by weight. 


STBIP 


Coolant inlet temperature, °F. 

(Note: STBIP is used only if inlet coolant tem- 


FOUL 


peratures are not read in or transferred.) 
Tube fouling resistance hr-f t 2 -°F/Btu 


ENHI 


Tube internal heat transfer enhancement factor. 



Card 4, (FORMAT: 5A1, 7F5.2, 415) 



GAS 


Type of noncondensible gas, must be one of AIR, 


SECFLG 


C02, or MIX. 

Bundle model indicator, gives number of 30° sec 
tors in one-half of a symmetrical bundle, or 
total sectors in a semicircular bundle. 


HFCDFL 


Symmetry indicator. 

If = 1, model will have symmetry about a ver- 
tical centerline 

If = 0, model will be a semicircular type. 


RADFLG 


Factor to increase bundle radius beyond that 
normally calculated by the code (will enlarge 
center void) . 
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BAFFLE 



FDAVE 



EXITFR 



0UTPUT 



INSTM 



ITRAN 



IFL0AT 



Flag to indicate presence of baffles at 2 
and 4 o'clock. 

Tube spacing parameter. Used in subroutine 
HETTRN to modify effect of condensate flood- 
ing . 

Target value for vented steam as a percent of 
input. If EXITFR = 0.0, then program will 
make a single pass and return the results 
without any adjustments. 

Control flag to restrict amount of detail 
printed out. 

0UTPUT =1.0 gives full detail of all sectors. 
0UTPUT = 0.0 prints only summary. 

Flag to direct program flow when converging on 
exit fraction (EXITFR) target. 

INSTM = 1, inlet steam adjusted. 

INSTM = 0, tube length adjusted. 

Flag to direct code to use outlet coolant tem- 
peratures from immediate preceding run (used 
only when multiple stages are run together) 
for inlet temperatures of current stage. 

ITRAN = 1, use transferred temperatures 
ITRAN = 0 , do not use transferred temperatures. 
Flag to give option of having detail output in 
floating point form; one gives floating point, 
zero gives fixed point. 
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IPNCH Outlet coolant temperatures will be punched 
for input if IPNCH = 0. 

Card 5, (FORMAT: 8F10.3) 

WSI Steam flow to condenser, lb/hr. 

WNCI Noncondensible gas flow to condenser, lb/hr. 

STSATl Inlet steam temperature, °F 

ALST Tube length, ft. 

ENH0 Tube external thermal enhancement factor for 
heat transfer. 

ENHF Friction factor enhancement factor for use in 
calculation of pressure drop in the steam. 

Card 6, (FORMAT: 4F1Q.3) 

ENHFI Friction factor enhancement factor for cooling 
water pressure drop 

TUBESW Specific weight of tube material, 1 /ft.^ 

FC Coefficient for enhanced- tube friction factor 

formula of the form f = FC* (Re cw /10 **FEXP 

FEXP See FC 

Card 7, (FORMAT: 215) 

IVAP, Switches to control vapor shear and inundation 

I RAIN 

effects. If set less than one these effects 
are neglected. 

Output 

Subroutine OUTP1 handles the principal output from the 
run. First, an image of the input deck is reproduced followed 
by an expanded version of this information. There next follows 
a series of statements defining the progress of the iterations 
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to achieve the specified exit fraction. In the baseline case 
(Appendix C) the tube length is successively adjusted until 
convergence is achieved in the ninth iteration. There fol- 
lows a table of data describing the condenser, cooler, and 
overall performance. Next, average inlet and outlet coolant 
temperatures are calculated for the condenser and the cooler. 
From these averages overall log mean temperature differences, 
DTCND2, DTC0L2, and DLT0T2, are calculated for the condenser, 
cooler and overall, respectively, using inlet steam tempera- 
tures. For the appropriate values of heat removed, the overall 
heat transfer coefficients, UPCOND , UPCOOL, UPAVG, are found 
for the condenser, cooler, and entire bundle. In addition, 
during the run, a row-by-row average of the heat transfer co- 
efficients (UBARW, UBARWC , UAVGW, respectively) is made, and 
the corresponding log mean temperature differences (ADTCND, 
ADTCLR , ADTOA , respectively) are back calculated. 

The option exists to print only the summary of results 
(OUTPUT = 0) , or to have the summary plus two pages of detailed 
results for each sector (OUTPUT = 1) . In the event of an error 
condition existing when the output subroutine is entered, the 
above calculations will be omitted and only the detailed output 
will be printed. 

Appendix C contains the summary output as well as the de- 
tailed output. In addition, a listing of definitions (not 
normally included) is inserted to facilitate the analysis of 
the detailed output. At the bottom of the summary output 
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(preceding the detailed output) are seven additional quanti- 
ties. These are defined as follows: 



BUNL2P 


Tube length for two-pass bundle, ft 


BUNVOL 


3 

Bundle volume, including tubes and voids, ft 


BUNWT 


Bundle weight, lb 


RETUB 


Cooling water Reynolds number 


FRIFAC 


Cooling water friction factor 


HDLOSS 


Cooling water head loss, ft H 2 0 


PMPPWR 


Cooling water pumping power, hp 




RESULTS OF COMPUTER EXPERIMENTS 



The following calculations examine applications of the 
MORCON code to departures from the baseline design. Several 
computer experiments are described in order to demonstrate 
some of the effects of internal and external enhancement of 
tubes contained within bundles. In addition, these examples 
serve to illustrate some of the influences that tube-to-tube 
interactions can have upon the performance of tubes in bundles 

Effect of Internal and External Enhancement on Bundle 



Leng th 
Figures 


8 a and 8 b illustrate the predicted effect of 



internal and external enhancement in terms of the length of 
tubes necessary to maintain the baseline heat duty. Constant 
heat duty is obtained by iteratively varying the tube length 
to obtain a fixed exit fraction of 0.05% of the baseline inlet 
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steam flow. Figure 8 includes manufacturer's data for en- 
hancement factors of some commercially available augmented 
tubes and it is seen that combined internal and external en- 
hancement can be expected to lead to bundle length (and volume 
and weight) reductions on the order of 35%. Table III gives 
information describing the augmented tubes used in this study. 



Table III. Characteristics of Some Commercially-Available 

Augmented Tubes ^ 



Tube 


Groove 

Geometry 


Thermal Enhancement 
Factors 


Head 

Loss 

Factor 


Pitch 

Dia 


Depth 

Dia 


Internal 


External 


KL 


0.543 


0.025 


2.14 


1.05 


2.03 


KM 


0.543 


0.040 


2.51 


1.14 


3. 26 


YL1 


0.250 


0.013 


1. 44 


1.35 


1. 56 


YL2 


0.250 


0.013 


1.44 


1.94 


1. 47 


YM1 


0.125 


0.0215 


2.70 


1.38 


6.02 


YM2 


0.125 


0.0215 


2.70 


1.94 


5.69 



Notes : 

1. Tubes KL and KM are products of Wolverine Division, 
UOP Inc. [25]. Tubes YL 1 , YL 2 , YM1 , and YM2 are 
products of Yorkshire Imperial Metals, LTD [26], 

2. Figures shown are ratios of head loss values to 
the smooth- tube case for the baseline conditions. 
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Figure 8. Effects 'of internal and external enhancement on 
tube length; baseline condenser. 
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The smooth- tube baseline is, as is typical, controlled 
in thermal resistance by the tube side film coefficient. For 
this reason internal enhancement is relatively more effective 
(about 2:1) than external enhancement. Figure 8b is a cross- 
plot of Fig. 8a and indicates the combinations of internal 
and external enhancement necessary for a given length reduc- 
tion. Here again it is seen that internal enhancement is 
justifiable to large values relative to external enhancement 
or, by way of example, external enhancement beyond 1.5 may 
not be warranted unless the internal film coefficient is en- 
hanced by a factor of at least 2. 

Variations in Thermal Resistance through the Bundle 

The importance of considering the limiting thermal re- 
sistance is further illustrated in Figs. 9a and 9b. Figure 
9a illustrates the row-by-row variation in internal and ex- 
ternal resistances for a condenser with no external enhance- 
ment. Without internal enhancement, the thermal circuit is 
unbalanced (dominated by the internal resistance) throughout 
the condenser. External enhancement in this case would clearly 
be gilding the thermal lily. With an internal enhancement of 
2.0, however, the two film resistances are relatively balanced. 
Such would be the approximate case for the KL tube of Fig. 8. 

If external enhancement is to be applied, on the other 
hand, such as in the case of YM1 in Fig. 8, additional internal 



46 



0 



to 

c • 

o to 

•H < 1 > 

+> a 

ftJ c 
■H ftJ 

P -P 
as to 

> *H 

to 

2 (D 
0 P 



P 

I rH 

>1 05 

-Q B 
I P 

£ CD 

o ^ 
es -p 



< j \ 

CD 

p 

a 

tn 

•H 

fP 




n^a/do 2 U Jq fr _0l *U 




rng/j 0 OL ‘a ‘soueisissu imiuam 

6 7 “ 



47 



Row number Row number 



enhancement is to be recommended. Such a case is illustrated 
in Fig. 9b where, again, the relative balance is obtained 
between the two film resistances by using an inside enhance- 
ment of E. = 3.0 and an outside enhancement of E =1.5. It 
1 o 

should be noted that in Fig. 9b the film resistances' have 
been lowered below the thermal resistances due to fouling 
(R c ) and tube wall resistance (R^) (resistance due to the 
non-condensible gas film is negligible in this case) . There- 
fore, further improvements would be most noticeable from re- 
duction in and/or 

Figures 9a and 9b also illustrate the complicated nature 
of the shellside resistance. The dashed lines indicate the 
distribution of external resistance in sectors I, IV, and VI 
of the condenser, with sector I at the top of the bundle (center 
ray at 15°) and proceeding clockwise in 30° increments to 
sector VI. In sector I a procession radially inward from 
row 1 to row 30 is in the direction of gravity and each suc- 
cessive row is inundated by the row above - thus the continuous 
increase in R in sector I. In sector VI, however, each row 
radially inward is subject to less inundation so that a gen- 
eral improvement is realized from the point of view of con- 
densate inundation. In every sector, however, the vapor 
velocity is decreasing in a radially inward direction and 
the beneficial effect of vapor velocity is correspondingly 
less in regions toward the center of the bundle. In the top 
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sector (sector I) condensate inundation and vapor shear ef- 
fects change in the same direction with resistance increasing 
radially inward. In the bottom sector (sector VI) the inunda- 
tion and vapor shear change oppositely and there is a general 
balance. It will be noted that near the bundle exit (entrance 
to the air cooler section) in sector VI the decreases in vapor 
shear effects no longer overcome the changes due to reduced 
inundation and the resistance begins to decrease. Sector IV 
is also illustrated in these figures and is seen to represent 
something of a mean behavior among the sectors. The spread 
between sectors in the external resistance, R q , is seen to be 
reduced in Fig. 9b because of reduction in the relative effect 
of vapor shear with external enhancement. 

Influence of Vapor Shear and Condensate Inundation 

Figures 10a and 10b illustrate the potentially deceptive 
nature of the influence of vapor shear and condensate inunda- 
tion in condenser bundles. Figure 10a examines the relative 
influence of the two effects upon required tube length and 
their dependency on internal enhancement. In this case the 
favorable effect of vapor shear is roughly balanced by that 
of inundation and in terms of required bundle length they are 
by no means independent of each other. Without enhancement 
the two effects are essentially equivalent so that the HE I 
method, which neglects them both, appears to be adequate (for 
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rather serendipitous reasons). With enhancement, however, 
when row-by-row thermal performance must be considered, the 
neglect of either vapor shear or condensate inundation can 
lead to erroneous results le.g., the neglect of vapor shear 
would indicate no improvement with an internal enhancement of 
about 1.7). Figure 10b shows that when external enhancement 
is large, these effects are less important because of the 
dominance of thermal resistances other than those due to the 
shellside film. 

Enhancement in Bundles of Various Diameters 

The number of possible comparisons such as those shown 
here is virtually limitless. As a final example, it might 
be useful to illustrate some of the influences to be ex- 
pected with enhancement when the number of tubes in the 
bundle and hence the bundle diameter is allowed to vary. 
Figures lla-llf illustrate this effect for a range of in- 
ternal enhancement factors with no external enhancement (Figs, 
lla-llc) and with E q = 3.0 (Figs, lld-llf) . In each set of 
figures the changes in bundle length, volume, and weight are 
shown for variations of + 20% in the number of tubes in the 
bundle with the cooling water velocity held constant at 8.0 
ft/sec. The reader will draw his own conclusions but a few 
interesting points may be highlighted here. 
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Comparing Figures a and d, b and e, and c and f, it will 
be seen that bundle length, volume, and weight may all be 
reduced through enhancement (E^ = E q = 3.0) by factors of 
about 38% with the number of tubes held constant at the base- 
line value of 3292. 

An increase in the number of tubes will require an in- 
crease in the bundle diameter while allowing, at the same 
time, a decrease in the tube length necessary to sustain the 
given heat load. These two effects are in opposition as re- 
gards bundle volume and weight and in the case of no enhance- 
ment (Figs, lla-llc, E^ = 1.0) the decrease in length (about 
27%) is sufficient to allow decreases in both volume (20%) and 
weight (13%) . At high levels of internal enhancement, how- 
ever, the allowable length reduction due to additional tubes 
is somewhat less and although volume reduction is still ob- 
tained the trend in bundle weight decreases and is eventually 
reversed with E^ = 2.5 and E q = 1. For example, the length 
reduction accompanying a 20% increase in tubes is about 11%, 
with a volume reduction of 5.4% and a negligible change in 
weight. At even higher values of E^ it is seen that a 20% 
increase in tubes leads to an increase in bundle weight - an 
effect that is opposite to that for the situation without en- 
hancement. Reversals in trends, such as these, are often 
impossible to predict if steam-side interaction effects are 
not taken into account. 
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If, on the other hand, it is desirable to reduce the 
number of tubes (for space considerations or to reduce cool- 
ing water flow rate, for example) a 20% decrease in the 
number of tubes will give a 10% reduction in bundle diameter. 
Without enhancement, this would require increases in length, 
volume, and weight of 90%, 70%, and 26%, respectively. If 
enhancement (E^ = E q = 3.0) is utilized, however, the per- 
centage increases in length and volume are only 15% and 6% 
while the bundle weight is decreased by about 7%. Thus, en- 
hancement can allow condenser modifications for "packaging" 
purposes that might otherwise be impossible. 

In the examples cited in Fig. 11, the condenser cooling 
water velocity has been held constant at 8.0 ft/s. For 
variable bundle diameters, the cooling water flow rate must 
therefore change in proportion with the number of tubes. 

These variations and constraints are used here for illustra- 
tive purposes and they represent only a few of the cases 
under consideration. It should also be remarked that here 
we have made only minor mention of the changes in cooling 
water head loss and pumping power that accompany the many 
variations shown in Figs. 8-11. Obviously, the relative 
merits of adjustments such as these depend upon the premiums 
to be placed upon the costs and benefits. If increases in 
bundle diameter can be tolerated, for instance, there are 
significant payoffs to be had in reductions in length, volume, 
weight, and cooling water pumping power. 
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SUMMARY 



The authors wish to again remind the reader that these 
calculations are most trustworthy when viewed as bases for 
comparisons. The extremely complex nature of the flow and 
heat transfer in condenser bundles leaves much to be accom- 
plished through further research. It is hoped that, with 
these precautions, the utility of a code such as MORCON has 
been illustrated. It is expected that future design com- 
parisons such as these will be expedited by the development 
of a numerical optimization code, based upon MORCON, to al- 
leviate the serendipity required to optimize such complex 
phenomena as condenser bundles. 

In closing, it is to be emphasized that several other 
performance factors and design variables, not mentioned here, 
are under study. The various tube-side and shell-side pres- 
sure drop legacies associated with enhancement are calculated 
and these constitute some of the "costs" to be balanced 
against the benefits of enhancement. A paramount factor 
yet to be understood is the extent to which fouling offsets 
the benefits of enhancement. Other design variables, such 
as tube material, tube diameter and wall thickness, condenser 
steam inlet pressure, and cooling water velocity and flow 
rate are under consideration. A change to titanium tubes 
will, in itself, generate a reduction in bundle weight by 
a factor of approximately 3 because of reduced metal density 
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and tube wall thicknesses. With all design variables left 
to vary within reasonable constraints, it is expected that 
improvements in tube bundle weight and volume may be realized 
in excess of those cited here and with acceptable penalties 
in terms of pressure drops and pumping powers. 

CONCLUDING REMARKS 



On the basis of work reported here and elsewhere, it is 
apparent that the continued advancement of the technology 
associated with surface condensers will require improvements 
in the understanding and definition of the flow conditions on 
the two-phase multi-component side of the problem. Whereas 
the prediction of the performance of single tubes - with 
vapor-side conditions well defined - has seemed adequate 
in the past, this situation may be attributed to conservative 
design attitudes and to the fact that thermal resistances 
other than those on the vapor side have dominated the overall 
thermal circuit. Advanced condenser designs that seek optimal 
performance (minimum weight, maximum efficiency, etc.) will 
require a level of analytical sophistication that exceeds 
present-day standards. 

A case in fact is that in which enhanced tubes are con- 
templated for use in condenser bundles. The prediction of 
condenser performance based upon single-tube theory and ex- 
periment cannot be expected to be reliable. Overall condenser 
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performance may exceed or fall short of such predictions 
depending upon the selection of the tubes themselves as 
well as their geometry and location in the condenser bundle. 

Computer-based analyses can be used to advantage in 
the estimation of the tube-to-tube variation of flow condi- 
tions within a condenser bundle. In this report, a condenser 
computer code of moderate complexity has been described and 
its use has been demonstrated in several sample cases. The 
results of these computer experiments have served to illus- 
trate the importance of internal detail in the prediction of 
the performance of condensers utilizing enhanced tubes. In 
addition, these studies have shown that judicious use of en- 
hanced tubes can lead to considerable improvements in con- 
denser performance as well as the realization of heretofore 
unattainable design options. 
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Appendix B. Input Data Deck for Baseline Design 
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